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Abstract 

A simple modification of the definition of the S-matrix is proposed. It is expected 
that the divergences related to nonzero self-energies are considerably milder with 
the modified definition than with the usual one. This conjecture is verified in a few 
examples using perturbation theory. The proposed formula is written in terms of 
the total Hamiltonian operator and a free Hamiltonian operator and is therefore 
applicable in any case when these Hamiltonian operators are known. 



1 Introduction 



The S-matrix, introduced in [HE], is a central object in scattering theory and its calcu- 
lation is of basic importance both in quantum mechanics and in quantum field theory. 
There are certain well known problems that one encounters in scattering calculations; 
one of them is that of the divergences caused by nonvanishing self-energies, which often 
appear in quantum field theory. In the present paper we concern ourselves mainly with 
this problem. 

A general definition of the S-matrix (which will be reviewed briefiy in section 12.11) 
can be given in terms of the total Hamiltonian operator H of the interacting physical 
system and a free Hamiltonian operator Hq. In quantum mechanics this definition can 
usually be applied without difficulty. In quantum field theory, however, one encounters 
divergences in perturbation theory arising from disconnected vacuum-vacuum diagrams 
and from radiative corrections on the external lines, if one takes in a straightforward 
manner the quadratic part of H as i^o (see e.g. |3j chapter 4. for a detailed explanation). 
The disconnected vacuum-vacuum diagrams and the radiative corrections on the external 
lines are known to correspond to shifts of the vacuum energy and of the masses of the 
particles (these shifts are called vacuum self-energy and self-masses, respectively) caused 
by the interaction part of H . 

The mentioned divergences are distinct from the usual ultraviolet or infrared diver- 
gences and are present whenever the vacuum self-energy and the self-masses are nonzero. 
They are not specific to quantum field theory either; similar divergences appear in any 
theory with nonzero self-energies, i.e. in any theory in which the eigenvalues of H are 
shifted with respect to those of Hq. 

This problem of divergences is usually solved in the framework of renormalization 
theory. If Hq and H are chosen appropriately, then the disconnected vacuum- vacuum 
diagrams cancel out entirely and the radiative corrections on the external lines get replaced 
by factors called field strength renormalization constants. 

In the present paper we propose a simple modification of the standard definition of the 
in and out states and thus of the S-matrix. The modification that we propose consists 
in including certain phase factors in the definition of the in and out states. Such a 
modification is allowed by the fact that physical states correspond to rays rather than to 
vectors, i.e. the phase of a state vector is not determined by the physical state that it 
represents. 

We suggest that with our definition the complications related to self-energy corrections 
are considerably milder than with the usual definition, and that our definition allows more 
general pairs of Hq and H operators than the standard definition. These are the main 
virtues of our definition. 

In section 12.11 we recall briefiy the standard formalism that we propose to modify. In 
the subsequent section [2T2] we present our modified formalism, and in section [2l3] we discuss 
some of its features; in particular its relation to a well known formula proposed by Gell- 
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Mann and Low ^i5j. The formulas in section are quite general; they can be applied in 
a wide range of fields of physics, regardless of symmetries and other particular properties. 
In section [3] we discuss various examples of different nature from quantum mechanics and 
from quantum field theory. These examples are aimed to provide some illustration for 
our definition and to demonstrate the cancellation of the divergences related to nonzero 
self-energies. We also verify that the results yielded by our formalism agree with the 
results that can be obtained by means of the standard methods. For our calculations we 
use perturbation theory, which is described in appendix [XI In appendix [B] we outline a 
regularization procedure for handling distributions (i.e. generalized functions) that occur 
in the calculations. Finally, we mention that in accordance with the nature of the subject 
of this paper we decided to cite mostly books and only a few papers. 



2 The general formalism of scattering theory 

2.1 A brief review of the standard definition of the S-matrix 

The S-matrix elements can be defined as 

Swv = {w , out\v , in) (1) 

where \v,in) and \w,out) are suitable states called in and out states. 

The in and out states are defined in terms of two self-adjoint operators H and Hq, H 
being the total Hamiltonian operator describing the scattering problem and Hq a free or 
reference Hamiltonian operator. The definition involves eigenvectors of Hq and the time 
evolution operator 

i 



Ueit2,ti) = Texp 



t2 



Hi,, it) dt 



ti 



(2) 



where 
and 



h 

Hj,,{t) = e-^l*le^^°*if7e-^^°* (3) 



Hj = H- Ho, (4) 
and the T in ([2]) denotes the time ordering: 

T[Ai{ti)A2{t2) . . . An{tn)] = v4fc^ (tfci )^fc2 (^^2 ) • • • (tfc„ ) , where ki,k2,...,K are deter- 
mined by the condition t^^ > t/.^ > ■ ■ ■ > tk„. It should be noted that U,{t2,ti) 
is unitary and it has the properties Ue{t,t) = I (where I is the identity operator), 
Ue(t2,ti)^^ = Ue(ti,t2) and Ue{ts,t2)Ue{t2,ti) = Ue{ts,ti). Au adlabatlc switching is 
included in the above definition of Ue{t2,ti), which serves as a regularization for the def- 
inition of the limits ^2 ^ oo or ti — > — oo. e > is the parameter of the adiabatic 
switching; 

hmU,{t2,ti) = U{t2M) = e*^o*2e-*^(*^-*^)e-^^«*\ (5) 
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The eigenvalue for an eigenvector \v) of Hq will be denoted by E^: 

Ho\v) = E'Jv). (6) 

The in and out states corresponding to an eigenvector |t>) of Hq are defined as 

\v,in) = lim lim UJO, —T)\v) (7) 
\v,out) = lim lim U,{0,T)\v), (8) 

SO the S-matrix elements are 

{w,out\v,in) = lim lim {w\Ue(T, —T)\v) . (9) 

e— >0+ T— >oo 

It is required that \v,in) and \v,out) be eigenvectors of H. We note that in this section 
and throughout the paper multiple limits are understood to be evaluated one by one from 
right to left. It is always the e — > O"*" limit which is evaluated last. 

In quantum mechanics the above definitions can usually be applied without difficulty. 
In quantum field theory, however, one usually finds that the e — 0"^ limits in (I7j), (IHl) 
and ([9]) do not exist if one takes the quadratic part of H as Hq. As we mentioned in the 
introduction, this divergence can be associated with a shift of the vacuum energy and of 
the particle masses caused by the higher degree terms in H. A standard method to deal 
with this difficulty is to modify (renormalize) the pair iJo? H used to produce the in and 
out states. 

For more detail on scattering theory we refer the reader to [31 El [71 El E], [THl [HI [121 [131 
[m [IS]. We note that instead of the adiabatic switching, other regularization methods 
such as taking the abelian limit are also used in the literature (see e.g. [IHl [TT]). 



2.2 The modified formulas 

In this section we present the modifications that we propose in the definitions ([71), ([HD of 
the |u,in) and \v,out) states, and thereby in the formula ([9|) for the S-matrix elements. 
We do not change 

For the in and out states we propose the following definitions: 

\v,in) = lim lim n^,,(-T)f/,(0,-T)|t;) (10) 
\v,out) = lim lim JT)UJO,T)\v), (11) 

e^0+ T->oo 

where Ilv,e(T) is a complex number of absolute value 1, i.e. a phase factor, given by the 
formula 



n.,.(T) = 





U,iT,0)\v){v 


Ue{0,T)\v) 


{v 


U,{0,T)\v) 



(12) 
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Ue{t2,ti) is given by the same formula ([2]) as in section [^?T1 The S-matrix elements are 
then given by 

= {w,out\v,in) = lim \imIlw^^(T)*Ily^^{—T){w\Ue{T,—T)\v). (13) 

The states \v, in) and \v, out) should be eigenvectors of H with eigenvalues that we denote 
by E^^in and E^^out- We do not require that = E^ and E^^out = E^. The eigenvalues 
Ey^in and E^^out can be obtained from the equations 



{v\HU,{0, 


-T)\v) 




~T)\v) 



, /n ' \ = ^^,in (14) 



{v\HU,{0,T)\v) 
hm hm Virrrn \ = " 

The definitions fllOj) and fllip are understood in the sense that if one wants to calculate 
a matrix element, for example, then one should do the calculation using 
limy^oo n^,^e(— 7')f/e(0, — T)|t>) and limT^Qo^v,eiT)Ue{0,T)\v) with e 7^ first, and take 
the e ^ 0+ limit afterwards. 

Often one has an orthonormal basis consisting of eigenvectors of Hq, where i 

labels the elements of the basis, and one is interested in the in and out states, S-matrix 
elements and energies corresponding to these basis vectors. We use the notation E^, Ei^m, 
Ei^out, ^i,e{T), Sij, etc. for this case. 



2.3 Discussion of the proposed formulas 

In this section we discuss the formulas proposed in section [221 In order to clearly identify 
the various remarks, we present them in the form of a numbered list. 

1. ) We stress that we do not require that the spectrum of Hq should be the same, 
even partially, as the spectrum of H. In particular, it is not required that the particle 
masses corresponding to Hq should be the same as the particle masses corresponding to 
H. Energy corrections, in particular the shift of the vacuum energy and of the particle 
masses can be calculated using (fT^ and (fTSl) . We expect that generally E^^in = E^^out- 

We also stress that in general neither \imT-,oo ^v,t{—T) nor lim-r^oo t^e(0, — T)|t>) in 
(ITOl) is convergent in itself in the e ^ 0^ limit (and a similar statement can be made 
about ffTTj) ). 

It is worth noting that the formulation in section [212] does not make use of symmetries, 
fields and the existence of a vacuum state. The formulas themselves are meaningful even 
for operators Hq, H which do not describe scattering. 

2. ) The phase factors Yi^^f^iT) and n^, e(— T) cancel out in the square of the absolute value 
of the S-matrix element (w, out\v, in): 

S^vS^j^ = {v,in\w,out){w,out\v,in) = lim lim {v\U^{—T,T)\w) {w\U^{T, —T)\v). (16) 

e^0+ T^oo 



5 



3.) The formulas ffTUj) . ffTTl) and ffT^ can be rewritten as 



, [7,(0, -T) It;) 
= ^lim jim ^X.,(-T) ^^^^^^^^_^^^^^ (17) 



jujO-ut) = lim lim \ Xy^^{T) 



U,{Q,T)\v) 



and 



'S'lui, = {w,out\v,in) = lim lim \ Xyj^^{T)\ Xy^^{—T) 



o+T^ooV ^'^^ '(t;|f/,(0,r)|t;) 

(w|f/,(T, -T)|t;) 



.o+T^ooV -'^^ "'^^ ^(u;|f/,(T,0)|u;)(i;|f/e(0,-T)|t;)' 

(19) 

where 

X,,,(T) = (t;|f/,(T,0)|i;)(i;|f/,(0,T)|t;). (20) 

The right hand sides of ffTUj) and ffTT]) written in the above form f ll7l) , fllSp are very similar 
to the formula introduced by Gell-Mann and Low in [4j (see also [5j), which serves to 
produce the vacuum state of H from the vacuum state of Hq. 

The right hand sides of ( |T71) and ( |T8|) and Gell-Mann and Low's formula (as written 
in [5]) differ in the factors Xy^^{—T) and ^/X^^^JT). These factors are real and are 
needed in order to maintain the correct normalization of the in and out states, which is 
important for the unitarity of the S-matrix. Of course, the formulas (fTTll and (fTSl) also 
differ from Gell-Mann and Low's formula in the type and range of eigenvectors of Hq to 
which they are applied. 

A main feature of Gell-Mann and Low's formula is that the phase of the vector that it 
produces from a vector \v) is fixed with respect to \v). This is true for our formulas ffTTl) 
and ffTSj) as well, since the factors ^/X^^J^—T) and ^JX^^^{T) are real. Let us introduce 
|fe(— T),m) and \v^{T),out) as 

\v,{-T),in) = n,,,(-T)?7,(0,-T)|t;) (21) 
\v,{T),out) = n,,,(T)f/,(0,T)K;). (22) 

We have 



X,,,(-T) = {v\v,{-T),in) (23) 



X,,,(T) = {v\v,{T),out). (24) 

Equation ( l23i) shows that if {v\v^{—T),in) ^ 0, then the overall phase of |fe(— T),m) is 
fixed in such a way that {v\v^{—T),in) is real. The same can be said about (12^ and 
|fe(T), out). 

The above feature of Gell-Mann and Low's formula is achieved by the denominator 
(f |[/e(0, — T)|w). If this denominator is not included, then the phase of the produced 
vector is generally not convergent as e — > (see |5j). 
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Gell-Mann and Low's formula is known to be applicable to discrete nondegenerate 
eigenstates other than the vacuum as well (see [H]), and it is not necessary for its appli- 
cability that the eigenstate of Hq to which it is applied have the same eigenvalue as the 
produced eigenstate of H. 

We can say that our proposal in this paper is to modify Gell-Mann and Low's formula 
(as written in [S]) by including the factors ^y Xy^^{—T) and ^yX^^^{T), and to apply it 
in this form to any eigenstate (not just to discrete nondegenerate eigenstates) with the 
purpose of producing general in and out states, which can contain any number of particles. 
As we mentioned above, the factors Xy^^(—T) and ^JJC^^^JT) are needed to maintain 
the correct normalization of the in and out states. The possible problems that can arise 
from the nondegeneracy of the eigenvector from which we intend to produce in and out 
states will be discussed briefly in section [3] in the fifth example. 

Concerning the requirement that the in and out states should be eigenvectors of H, 
Gell-Mann and Low's formula is known (see [HE]) to produce an eigenvector of H if the 
e ^ O"*" limit that it involves exists, and it is natural to expect that this remains true if 
the formula is applied to eigenvectors \v) of Hq which are not discrete and nondegenerate. 

4. ) Let us assume that \v) is a discrete nondegenerate eigenvector. In a perturbative 
framework \v,in) and \v,out) will also be discrete nondegenerate eigenvectors, and they 
will belong to the same one- dimensional eigenspace, i.e. \v,in) = a\v,out) with some 
complex number a. {v\v,in) and {v\v,out) are both positive real numbers (they can be 
assumed to be nonzero in a perturbative framework), hence it follows that \v, in) = \v, out) 
and Syy = {v\v) (see also [5]). This result holds in particular for the vacuum in quantum 
field theory. We expect that it also holds for one-particle states in Poincare symmetric 
theories. 

5. ) In the case of Poincare symmetric theories we also expect that for a multi-particle 
state \ki, ^2, kn) of Hq the total energy of \ki, /c2, kn, in) and \ki, ^2, kn, out) has 
the properties that Ek^^k2,...,kn,in = -£'fci,fc2,...,fc„,o«t and £^A;i,fc2v,fcn,m — -Sr2o,m = ^^=iiEki,in — 
-E'ncm), where Qq denotes the vacuum state. The latter additivity property is expected 
because \ki, k2, kn, in) and \ki, /c2, kn, out) should be multi-particle eigenstates of H. 

6. ) We expect that instead of the exponential switching function e~^'*' other functions 
can be used as well. For example, a linear switching function could also be used. The 
term adiabatic switching does not, of course, refer to any physically real switching of the 
interaction. In this paper we shall not discuss the precise role of the adiabatic switching 
prescription. A feature of this method that is worth noting is that the unitarity of the 
time evolution operator is preserved for all values of e. 

7. ) In quantum theory the vector representing a physical state is fixed up to a phase 
factor only (in other words, physical states correspond to rays), so in general we can say 
that the phases of the in and out states are not determined: one could multiply the right 
hand sides of (ITOl) and (|TTi) by any phase factors Hy^in and Tiv,out- 
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8.) Since the in and out states are eigenstates of H, 



S^, = (25) 

holds if Eui^out 7^ Ey^in- This property expresses the conservation of energy. 

9. ) In the case of some Hq and H it can happen that ( fTOl) and ( |TT1) do not produce 
eigenvectors of H for certain eigenvectors of Hq. One of the situations in which this can 
be expected to occur is when Hi makes some particles unstable. 

10. ) One encounters infrared and ultraviolet divergences in many theories. In this paper 
we do not discuss these divergences; they can be handled in the same way as usual. 

11. ) Our formalism allows more general pairs of Hq and H than the standard formalism 
(i.e. ([7]), (IH]), (EI)), nevertheless the right choice of if, Hq and {\vi)} for a given physical 
system should be subject to consideration in general. We do not discuss the physical 
interpretation of Hq and {|fi)} in this paper; we regard them as auxiliary quantities 
which are used to produce, by means of the formulas (ITOl) and (fTTl) . suitable in and out 
states. 

12. ) In certain problems, e.g. in rearrangement scattering, one can have two different free 
Hamiltonian operators Hq and if^f ; one for the incoming states and one for the outgoing 
states (see e.g. [H]). Even if = Hq = Hq, it might be necessary to take two different 
sets of basis vectors and to define the in and out states. For simphcity, we 
restrict ourselves to the case when Hq = H^ = Hq and a single set of eigenvectors is used 
to define both the in and the out states. The extension of the discussion to the general 
case is straightforward. 

13. ) By orthonormality of the eigenvectors \vi) of Hq we mean 

{vi\vj) = 6{i,j), (26) 
where S{i,j) is a Dirac-delta; S{i,j)=0 ii i ^ j and 

d3 5{t,3)g{j) = g{{), [ dt5{t,j)g{t) = g{j), (27) 



where g{i) is any test function. The integration J di over the index set is understood in a 
general sense; it may include summation over discrete parts of the spectrum, for example. 
In accordance with this, the Dirac-delta S{i,j) is also understood in a generalized sense. 
The name eigenvector is often reserved to proper eigenvectors with finite norm, and eigen- 
vectors belonging to continuous parts of the spectrum are called improper eigenvectors. 
In this paper we use the term eigenvector for both proper and improper eigenvectors. 

14.) {vi\vi) is often not a finite number and {vi\vj) has to be regarded as a distribution. 
When one applies perturbation theory or other methods to calculate (JT3l) and the other 
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quantities in section 12.21 this leads to expressions which are not well defined, in princi- 
ple, since ffT^ will involve products and quotients of distributions. One way to overcome 
this difficulty is to introduce a regularization of the eigenvectors \vi). Such a regulariza- 
tion prescription is outlined in appendix [Bl In some cases, however, it is sufficient to 
manipulate with {vi\vj) in a formal manner, and the regularization can be avoided. 

15.) As is well known, the nonzero S-matrix elements (and also the components of the in 
and out states) are usually not finite numbers, but the integrals / di dj (f){a,i)*(j){b, j)Sij, 
where 0(a, i) and 0(6, j) are suitable functions, are finite, i.e. the S-matrix is a distribution 
(and a similar statement can be made about the in and out states). The mentioned 
integral is equal to where |$a) and are wave packets of in and out states: 

= / di(f){a,i)\vi,out), = / dj (t>{b, j)\vj,in) . 



16. ) In section 12.21 the quantities called S-matrix elements are defined for a certain 
set of in and out states that correspond to eigenvectors of Hq. However, a matrix 
that could be called S-matrix was not defined. In order to be able to define a com- 
plete S-matrix we assume that an orthonormal set of vectors {l^j)} is chosen; then the 
matrix constituted by the S-matrix elements Sij can be called the S-matrix. The S- 
matrix elements for arbitrary superpositions of the in and out states \vi,in) and \vi,out) 
are determined by linearity: let |$a) and be superpositions of in and out states: 
l^a) = / di(f){a,i)\vi,out), = / dj (f){b, j)\vj,in) . The S-matrix element Sab is then 
Sab = {^a\^b) = / di dj (f){a, i)* (f){b, j) Sij . It is important to note here that the \v) ap- 
pearing in (fTUl) and (ITTI) should always be an eigenvector of Hq, and that (|TU|) and (fTTI) 
are not linear in \v). 

17. ) It is easy to verify that the unitarity of the S-matrix, which is expressed by the 
equation f dk Sik{Sjk)* = {vi\vj), is not affected by the phase factors included in (fTUI) and 



In this section we discuss some results that we obtained in certain specific models by means 
of the application of (|T0l) -( fT5l) . H takes the form H = Hk + fi'-f^mt in all cases, where g 
is a coupling constant. We studied five examples, which are the following. The first one 
is the general case when has a discrete nondegenerate and finite spectrum. Although 
the Hamiltonian operators of this type do not describe scattering, from a technical point 
of view we found it interesting to consider this case. The second example is the scattering 
of a single particle in short range spherical potentials in three dimensional space in the 
framework of nonrelativistic quantum mechanics. Scattering in a Dirac-delta potential in 
one dimensional space is also briefiy discussed. The third example is the scattering of a 
massive relativistic particle on a defect in 1 + 1 spacetime dimensions in the framework 
of quantum field theory. The particle is a real scalar boson and the defect is localized at 
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a; = 0. The free particle Hamiltonian operator is 



1 /"°° 

Hk = - dx [: {dt^f + {d^^f + m2$2 :]; (28) 

J -oo 

the interaction term is 

i^int =: $(0, 0)' : . (29) 

The fourth example is the theory in 3 + 1 spacetime dimensions. The free particle 
Hamiltonian operator is 

HK = \j d^x [: {dt^f + {d^^f + m2$2 :]; (30) 



iJint = / d=^x : <l>^ : . (31) 



the interaction term is 



The fifth example is the case when Hk has a discrete, finite, but not necessarily nonde- 
generate spectrum. 

We took Hk (i.e. the constant part of iif as a linear function of g) as the reference 
Hamiltonian operator Hq in all of the examples. 

We used for our study the framework of the ordinary perturbation theory, as described 
in appendix |Al In the first and third cases we did calculations up to second order in g\ 
in the fourth example we also considered the third order. We considered all orders in 
the second example. Feynman diagrams can be associated with the various terms also in 
ordinary perturbation theory, if it is applied to quantum field theories like the theory. 

In the first example we found that our formalism reproduces the results that can be 
obtained by Rayleigh-Schrodinger perturbation theory (see e.g. p!6]) and gives the result 
that the S-matrix is the identity operator, in accordance with the remarks in 4. in section 

As the second example we discuss potential scattering. We expect that the modified 
formula fll3p and the standard formula ^ yield the same S-matrix in a large class of 
potential scattering problems. In order to support this conjecture we will now discuss 
the scattering of a single particle in short range spherical potentials in three dimensional 
space. We consider all orders of perturbation theory. 

The basis vectors of the free particle Hamiltonian operator = are denoted 

by |k); their wave function is 

Ik) = ^ e'*^" (32) 



The energy of |k) is k 



2. 



i/x|k) =k2|k). (33) 
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The interaction Hamiltonian operator is 



^^int = ^(x), (34) 

where ^(x) is the potential in which the scattering takes place. We assume that l^(x) is 
spherically symmetric, continuous and there exists a number p > such that it falls off 
at least as fast as |x|~^~'' at infinity (i.e. |V"(x)| < c|x|~^~'' if |x| > ro, where c and tq are 
suitable constants). Certain singularities at x = and at other points could be allowed 
as well (see chapter 2 of |lDj). The matrix elements of ifint are (see ( 153|) for the notation 
on the left hand side) 

<^-''^> - (2^/^^^ e*.-.>"nx) -^,fr',r ^^Vir), (35) 

where r = |x| and q = |ki — k2|. 

In the following we will show that limT^oo nk,e(— T) = 1 and limT^oo nk,e(T) = 1, 
which imply that ( fT3i) and ([9]) yield the same S-matrix. We consider limr^oo nk,e(— T) 
explicitly; the case of limT^oo nk,e(T) is very similar. We use the notation Ckik2,K,- 
the coefficients in the Taylor series for limT^oo(ki|?7e(0, — T)|k2) (see ([SI])), e is kept fixed 
at some arbitrary positive value. 

Under the above assumptions for ^(x) the matrix element (kik2) is finite for any ki 
and k2 and |(kik2)| is also bounded. As one can see from ( |35l) . |(kik2)| falls off for large 
q at least as fast as 1/q. The coefficients Ckik2,K,-) ^ = 1)2, ... , given by the integrals 
f l3U|) . are also finite for any finite value of e if one introduces a momentum cutoff. Power 
counting shows that the integrals giving these coefficients are finite without a cutoff as 
well, i.e. they are not ultraviolet divergent. 

We introduce the regularized vectors (see appendix [B]) as 

|k,/i) = j d^k' x^iikl,K)x^iik2,k'2)x^iih,K)\^')y (36) 

where ki, k2, are the components of k, /x > is a small real number and x^i{k, k') is 
the function given by 

X^.{k,k')= 2^ ii \k-k'\<k,^^ (37) 
Xf,{k,k')= if \k-k'\>koH, (38) 

where fco > is an arbitrary fixed constant. We use the notation Ckik2,K,M,- for the 
coefficients in the series for limT^oo(ki, /i|t/e(0, — T)|k2,/i). We have 

Ckik2,K,M- = / d^k'id^ka Xf,i^i,K)Xf,(M2,K)C]^[y:,,t^- , (39) 
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where x^l{K'^') = XtJi{ki,K)xi,{k2,k'2)Xii{h,k'^)- A coefficient Ckika^./^ where k > 0, 
has the hmit Ckik2A,- as yU ^ 0, which is finite. From now on we will take ki = k2 = k, 
since this is the relevant case for nk,e(— T). The zero order term in the Taylor series for 
limT^oo(k,/i|f/,(0, -T)|k,/i) is 

which goes to infinity as — 0. The Taylor series for limr-^oo nk,e,^(— T) can be ob- 
tained from that for limT^oo(k, /i|t/e(0, — T)|k, /i), taking into consideration (IH^ . The 
zero order term in the Taylor series for limT^^oo ^'k,e,ii{~T) is clearly 1. The coeffi- 
cients of the further terms in the series will be polynomials (with zero constant term) 
of Ckk,K,/i,-/ (k, //|k, /i) and C*j^k,K,/i,-/ (k, /^|k, yu), where k can take positive integer values. 
However, Ckk,K,M (k, yu|k, /i) and CjJk,^,/^ -/ (k, A^lk, /i) go to zero as /i ^ 0, since Ckk,K,At,— 
has finite limit and (k, /i|k, /i) goes to infinity. Thus all coefficients in the Taylor series 
for limT^oo nk,e,^(— T) go to zero as /i — > 0, only the zero order term remains, which is 
1. This completes our derivation of limy^oo nk,e(— T) = 1. In summary we can say that 
the essential point in the derivation is that Ckk,K,Af,-/ (k, /i|k, /i), k = 1,2,..., go to zero 
as /i — 0. 

A similar argument as above can be applied to show that the formulas (fT4l) and (ITSl) 
give zero for the corrections to the energy eigenvalues. We discuss (fT^ only; the case of 
(fT3|) is very similar. Using H = + gH^^x. we have 

^. (k,/x|g^,(0,-T)|k,/i) ^ (k,/x|gK^,(0,-T)|k,/i) 
T™ (k,/x|?7,(0,-T)|k,/x) T™ (k,/i|f/,(0,-r)|k,/x) 

+ ,hm (k,/^|i^.n.^e(0,-T)|k,;.) 

+ ^T^^ (k,Hf/e(0,-T)|k,/i) ^ ^ 

The limit of the first term on the right hand side as /i — > is the unperturbed energy 
i?^ = k^. In the second term the numerator has a Taylor series in which all the coefficients 
have finite limit as /x — >■ 0. The coefficients in the Taylor series for the denominator also 
have finite limit as /i ^ with the exception of the zero order term (k, /i|k, /i), which 
goes to infinity as l//i^. All the coefficients of the Taylor series for the second term will 
therefore go to zero as /i — 0, thus the second term is zero, and -Ek.m = -E'k = k^. 

Very similar derivations to those presented above can be applied also, for instance, 
for the scattering of a particle in the Dirac-delta potential in one dimensional space. In 
this case the Hk operator is the free particle Hamiltonian operator Hk = the 
interaction Hamiltonian operator is ifjnt = 5{x). We note that the S-matrix takes the 
form 

Sk,M = Kki - k2)T{k2) + 6{ki + k2)R{k2), (42) 
and the exact expressions for T{k2) and R{k2) can also be found; they are 

T{k2) = ,, 'j^'' , , R{k2) = ,, ^J'^ , . (43) 
^ i\k2\+g/2' ^ ^' i\k2\+g/2 ^ ' 
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We also note that in perturbation theory there are infrared divergences in this model for 
in and out states with zero momentum. 

Long range potentials, such as the Coulomb potential, present similar problems for 
the application of both ( fT3l) and the standard formula ([9]). In the case of the Coulomb 
potential, for instance, it is not immediately obvious what the matrix elements (kik2) are, 
since the integral on the right hand side of (1351) is not convergent for large r. A possible 
method of handling long range potentials is to introduce a shielding, i.e. to approximate 
them by short range potentials, as described in detail in [10]. This method is suitable for 
(fT3l) as well. 

In the third example the interaction term breaks the Poincare symmetry of the free 
particle, so the 1-particle S-matrix 5*^1,^2 is nontrivial. Sfc^^^j is given by the same ex- 
pressions f l42l) . fHHj) as the S-matrix for the scattering in the delta potential mentioned 
above, and our perturbative results were in agreement with this exact expression. For the 
vacuum- vacuum S-matrix element we obtained Sn^^Ug = l + 0{g^), in agreement with 4. in 
section [231 Concerning the energy levels, we found that all states, including the vacuum, 
get the same finite correction, the first order part of which is zero. We considered parti- 
cles in in and out states with nonzero momentum only, since at zero momentum infrared 
divergences occur. We note that this model was also studied in a different approach in 

Turning to the fourth example, we assumed, in accordance with the perturbative 
framework, that the theory describes an interacting massive relativistic boson. It is 
well known that in this model one encounters ultraviolet divergent loop integrals in per- 
turbation theory. Such divergences do not have significance for the present study and 
they do not cause much difficulty, therefore we were not concerned with them. One can 
introduce a momentum cutoff if one wants to have a regularization of these divergences. 
Radiative corrections on external lines appear first at third order, therefore we calculated 
S-matrix elements up to third order in this example. The results can be compared with 
those that can be obtained by means of the usual bare perturbation theory (in which H is 
not renormalized) . There are some differences between the presentations of standard per- 
turbation theory in quantum field theory in various textbooks; although these differences 
are not essential in principle, we note that we used [3j as reference. 

For the vacuum- vacuum S-matrix element we obtained Sq^^^q^ = 1 + O^g"^), in agree- 
ment with 4. in section [^751 For the one-particle S-matrix we got 5'ki,k2 = '^'^(ki — k2) + 
0{g^), again in agreement with the expectation stated in 4. in section As regards the 
two-particle S-matrix elements >S'k3,k4;ki,k2 5 we obtained a result which is identical to the 
result that one can obtain in usual bare perturbation theory (see e.g. [3] for a description 
of bare perturbation theory). The result for 5'k3,k4;ki,k2 is a formula that contains some 
loop integrals, which are ultraviolet divergent, of course. 

We also calculated the correction to the mass of the particle to second order. This 
correction can be obtained from the energy levels in the following way: The mass m of 
a one-particle state is given by m = \JE'^ — k^, where E = E]^ — Eq^ is its energy and 
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(a) 



(b) 



Figure 1: Second order diagrams for the mass correction. 



k is its momentum. iJ^ and E^g denote the eigenvalue of the one-particle state and the 
eigenvalue of the vacuum state with respect to H. We denote the expansion coefficients 
of E and m as follows: 

E = E^ + gE^ + g^E'^ + ... , m = mo + gmi + g^m2 + ... . (44) 

Similar notation as for E applies to the expansion coefficients of E^ and Eq^ as well. The 
series expansion of m in terms of these coefficients up to second order in g is 

m = mo + g + ... , 45 

where mo = a/ {E^Y — k^. It is important in ( l45l) that mo 7^ 0. In the present example 
-E^ = 0, therefore (H5l) takes the form m = mo + g^E^E'^ /mo. E^. and E^^, and thus 
E"^ = E^ — Ef^^, can be calculated using ( |78l) . The terms contributing to E"^ are associated 
with the graphs shown in figure [TJ The contributions corresponding to the graph [Hb are 
equal to Ef^^ and so are cancelled out entirely by the subtraction of Ef^^. The terms 
contributing to E'^ will therefore be those corresponding to the graph [Ha. The result 
that we obtained for m2 agrees with the result that can be obtained in bare perturbation 
theory in the standard formalism, in which the particle mass is given by the location of 
the pole of the propagator. The result for the mass correction is a formula that contains 
an integral which is, of course, ultraviolet divergent. 

Concerning multi-particle states, it is not difficult to verify that the additivity property 
mentioned in 5.) in section [231 holds. Only those graphs give contributions to the energy 
of a multi-particle state which are totally disconnected, i.e. which do not describe any 
interaction between particles. The contributions of those graphs which describe some 
interaction other than self-interaction turn out to be zero in the /i — > limit (where /j, is 
the regularization parameter for the eigenvectors of Ho, see appendix [B|) . 

The fifth example is included in order to illustrate certain technical points related 
to the degeneracies of the eigenvalues of Ho. Let us assume that \vj) is a degenerate 
eigenvector. Applying perturbation theory to calculate {vi\vj,in), at first order we find 
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that the e ^ 0+ hmit is convergent only if (f j|ifint|fj) = for all values oil ^ j for which 
= Ej. This condition is known from the degenerate Rayleigh-Schrodinger perturbation 
theory and can be satisfied by choosing the basis vectors within the eigenspace of Hq 
containing \vj) in such a way that the matrix {vi\Hiat\vj) be diagonal within this eigenspace 
of Hq. Assuming that this condition is satisfied, at second order we find that the e — 0"*" 
limit exists if and only if 



E 



{Vi I if int I Vm ) {Vm \ ^^int 1^ - Q (^45 ) 



E^ -^0 

m J 



holds for all values oi i ^ j for which E^ = Ej, where the summation has to be done 
over those values of m for which E^ ^ Ej. This condition is not exactly the same as the 
one in Rayleigh-Schrodinger perturbation theory; in the latter case, the equation (H6l) is 
required to hold only for those values oi i ^ j for which E^ = Ej and Ef = Ej {E} and 
Ej denote the coefficient of the first order corrections to E^ and Ej), i.e. for which \vi, in) 
remains degenerate with \vj,in) at first order. While the weaker condition of Rayleigh- 
Schrodinger perturbation theory can always be satisfied by choosing suitable basis vectors, 
the stronger condition that we have found cannot always be satisfied if the degeneracy 
of the eigenspace containing \vj) is broken at first order. This problem can be handled 
by shifting those eigenvalues E^ of Hq for which E^ = Ej but E} ^ Ej, i.e. by taking 
Ho = Hk + where SEi are the energy shifts. There is some freedom in the 

choice of 6Ei] one can take 6Ei = g{El — Ej), for instance. Hj = H — Hq also changes 
as a result of changing Hq] this should be taken into consideration in the perturbation 
theory. We expect that at higher orders we would find conditions similar to P6|) . which 
could also be handled in the same way as fHUl) . In quantum field theory or in quantum 
mechanics one could also find conditions which are analogous to those described above. 
These conditions, which we do not discuss in detail, are often milder for those eigenvectors 
of Hq which belong to continuous parts of the spectrum. If, for instance, {vi\Hirit\vj) = 
were necessary in any case when i ^ j and Ef = Ej, then nontrivial scattering would 
be impossible at first order. In some cases, for example when an interaction breaks the 
mass degeneracy of a multiplet, a suitable modification of an initial choice of basis vectors 
may be necessary; in some complicated quantum field theories a modification of Hq may 
be necessary as well. We can say that with our definition of the in and out states the 
divergences of the e — *• O"*" limit are related to certain kinds of breaking of degeneracies of 
the eigenvalues of Hq by the interaction. 

The calculations needed to obtain the above results are largely straightforward, nev- 
ertheless lengthy in some cases, especially in the $^ theory. The precise treatment of the 
e — s> 0+ limit is essential; in particular the coefficients of e in the denominators in ( l56l) . 
( I59l) . ( l63l) must not be changed (see appendix |A] for some further comments). In the 
second, third and fourth examples it is also important to keep in mind that the S-matrix 
elements and the components of the in and out states are distributions in general. One 
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also encounters the problem mentioned in 14. in section I2.3[ As we said there, one can 
introduce a regularization to handle this problem. We did this in the second example, 
otherwise it was satisfactory in our calculations to handle the Dirac-deltas in a formal 
manner. 

In the first, third and fourth examples the {w\Ue(T, —T)\v) , Ilw,e{T)* and Ilv,t{—T) 
parts of (fT^ are not individually convergent in the e ^ 0"*" limit (after T — > oo, of 
course). Up to the orders that we considered we found that these parts contain terms 
which diverge as 1/e in the third and fourth examples; in the first example terms di- 
verging as 1/e^ also occur. In the third and fourth examples these divergent terms in 
{w\Ue{T, —T)\v) can be associated with diagrams having vacuum- vacuum diagram parts 
and radiative corrections on the external lines; the latter occur only in the fourth ex- 
ample. The phase factors Ilw,e{T)* and Ily^^{—T), however, give rise to further terms 
in the product Ilw,eiT)*Ily^^{—T){w\U^(T, —T)\v), which turn out to have the effect that 
the divergences that would occur in the e — 0"^ limit are eliminated; all diagrams with 
vacuum-vacuum parts cancel out and the radiative corrections on the external lines are 
replaced with a constant factor, which turns out to be equal to the square root of the 
field strength renormalization constant of the standard LSZ formalism. This constant, 
which appears in the fourth example only, is given by a formula which contains ultraviolet 
divergent integrals. The reader is invited to carry out these calculations to see how the 
cancellations and the replacement take place. In summary, we arrive effectively at the 
usual rules of perturbation theory, without the need to modify the choice Hq = Hk- As 
we mentioned in the introduction and in section [^TTj if one wanted to apply ([H]), then 
one would have to modify Hq (assuming that H is not changed) because of the vacuum- 
vacuum parts and the radiative corrections on the external lines. We also verified in all 
of the examples up to second order that the e — O"*" limit of the components of the in 
and out states (i.e. the scalar products of the in and out states with the free states) is 
convergent up to second order (the ultraviolet divergences in the $^ theory are present, 
of course, for any values of e). The eigenvalue equations for the in and out states are 
also satisfied up to second order. Since the eigenvalue equation is an equation of vectors, 
we took the components of these vectors with respect to the free states and verified the 
equality of these components. As we mentioned earlier, it is important in the case of the 
second, third and fourth examples to keep in mind that one deals with distributions. 

Finally, it is also worth emphasizing that the results 5'f7o,no = 1 ajid, in the theory, 
S'ki,k2 = "^^(ki — k2) could be derived in our formalism in a straightforward way; it is not 
necessary to assume or to require that these equations hold. 

4 Conclusions 

In this paper we proposed a modification of the traditional general definition of the in 
and out states and of the S-matrix elements. The main components of our proposal are 
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the following: phase factors are included in the definition of the in and out states in order 
to prevent oscillations occurring due to nonzero self-energies; an adiabatic switching is 
used; and the requirement that an in or out state has to have the same energy as the 
corresponding free state is dropped. The divergences related to nonzero self-energies 
are considerably milder with the modified definition than with the usual one. With the 
modified definition divergences can occur due to certain breakings of degeneracies by the 
interaction. 

Using perturbation theory, we verified in certain quantum mechanical and quantum 
field theoretical models that our definition yields the same S-matrix elements as the usual 
definitions. In the quantum field theoretical examples we did the verification up to second 
and third orders. For those cases of potential scattering that we discussed we presented 
a derivation of this result to all orders. In the quantum field theoretical examples it 
was not necessary to do any renormalization to remove the divergences related to nonzero 
vacuum self-energies and self-masses. We also calculated energy and mass corrections and 
obtained results which are in agreement with expectations and with the results that can 
be obtained using standard definitions. We expect that these results extend to arbitrary 
orders of perturbation theory, and similar results hold for other models in quantum field 
theory as well. 
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A Perturbation theory 

In this appendix we describe ordinary perturbation theory, which is our tool in section 
El We include this part in the interest of completeness and because the coefficients of e 
in the denominators in ([55]), (IMD, (EHD, (EH]), (JM]), ([63]), (JM]), (JET]) and (EH]) differ from 
those that can be found in similar formulas in a large part of the literature. However, 
these coefficients should generally be treated precisely in the calculations; otherwise one 
obtains wrong results in some cases. The coefficients of e are characteristic of the adiabatic 
switching prescription and of the exponential switching function. 
We assume that H takes the form 

H = Ho + gHint, (47) 

where g is a coupling constant and Ho is the reference Hamiltonian operator. The Hj 
operator is thus Hj = gH^at- We also assume that an orthonormal basis of eigen- 

vectors of Hq is given. The perturbation series which we shall deal with will be Taylor 
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series in g. We assume that Hq, Hi^t and is independent of g. If these quantities 

depend on g, then the series written below will not be Taylor series in g; in order to obtain 
the power series one should expand the various terms of these series into powers of g and 
collect the terms proportional to the same power of g. 

A series expansion for {vi\Ue{t2,ti)\vj) can be obtained in the following way. One 
expands U^(t2,ti) into the Dyson series 



Ue{t2,ti) = Texp 



t2 



gH,^Ut) dt 



N=l 

oo 



Hint,e{TN)] 



'tl 

^ ft2 PTl /•T2 

'^Jls^ii) / dr,/ drj dr,... 

j^^l \ / Jti Jti Jti 

f T]V_1 

dr7vifint,,(n)^^mt,e(T-2) 

'tl 



HinU^N), (48) 

where Hi^^^t) = e~^'*'e^^°*ifinte~^^°*. One inserts unit operators in the form 



dm \Vm){Vr, 



(49) 



between the -f/'int,e(T)-s in 

{Vi\Hi,,t,e{^l)Hint,ei^2) ■ ■ ■ HintA^N)\Vj) = 

dmidm2 . . . dm7V_i {Vi\Hint,e{Tl)\Vmi){Vmi\Hint,e{'^2)\Vm2) ■ ■ ■ (^^m^r_l | -f^int,e (^iv) ) , 



and then performs the integrals over the r-s, which is possible since 



{Vn\H;^t,e{T)\Vm) = (t;„ | -f^int | f m) exp 



(50) 



(51) 



is a simple exponential function. After the integrations over the r-s are performed, the 
integrations J dmidm2 . . . dm7v-i over the intermediate states still remain in the formula. 

Taylor series that we obtained for \imT^oo{vi\Ue{0, —T)\vj), \imT^oo{vi\Ue(T,0)\vj) 
and \imT^oo{vi\U^(T, —T)\vj) are presented below. These formulas arise after performing 
the integrations over the r-s. 

The series written below for {vi\U^{ti,t2)\vj) , \imT^oo{vi\U^{0, —T)\vj), 
li'mT^oo{vi\U6iT,0)\vj) and limT^ooivilUeiT, —T)\vj) can be used to obtain the power 
series for the various quotients and products in section [221 We list the coefficients of the 
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series that we obtained for the in and out states, for the S-matrix, and for the eigenvalues 
of H up to second order. Coefficients for some further quantities are also given. 

We remark that the matrix elements {vi,in\A\vj,in) , {vi,out\A\vj,out), 
{vi, in\A\vj, out) and out\A\vj, in) of any operator A between in and out states can also 
be calculated perturbatively if the matrix elements (^11741^^) of A between the eigenstates 
of Hq are known. For instance, 

{vi,in\A\vj,in) ^ J dmidm2 (t-i, m|f^i)(t'„JA|t'„2)(t'„2|T;j, m), (52) 

and here the perturbative results for {vi,in\vjni) and (vm^lvjjin) can be used to obtain a 
power series for {vi,in\A\vj,in). 

We begin our list of formulas with formulas for \imT-^oo{vi\Ue{0, —T)\vj), 
limT-*oo{vi\Ue{T,0)\vj) and limT-^oo{vi\Ue{T,—T)\vj). We use the notation 

P{ab) = '-{El - El), {ah) = {a\H,,,\b). (53) 
A.) For \imT-*^{vi\Ue{0, —T)\vj) we obtain 

]imJv,\U,{0, -T)\v,) = E ( X ) ' ^^^^ 



where 



Cijfl- - {Vi\vj), Cij^i- - -^^rry^, (55) 



a 



I 



, , , {imk^i){mk-imk-2) ■ ■ ■ {ni2nii){'mij) 

mi m2 . . . ruk-i ^^^^^^ ^ ke][P{mk-i3) + [k - l)e] . . . [P(m2j) + 2e][P{m,j) + e] ' 

(56) 

B.) For limT^(x,(i^i|^e(^, 0)|fj) we obtain 



OO / . \ fc 

^lim (^;,|C/,(r, 0)1^;,) = E ( ^ ) 3'C,^,k,+ , (57) 

k=0 ^ ^ 



where 



{ij) 

Cij,Q,+ = {vi\vj), Cij,i,+ = p^-^yip^' (58) 



J 



{imi){mim2) . . . {mk-2mk-i){mk-ij) 

dmidm2 . . . dmjk_i 



[P(mii) + e][P{m2i) + 2e] . . . [P{mk-ii) + {k - l)e\[P{ji) + ke] ' 

(59) 
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C.) For limT^^{vi\Uf:{T, —T)\vj) we obtain 

^^Jv^T, -T)\v,) = E ( Y ) ' 

h—n \ / 



fe=0 



where 



Dij^o = {vi\vj) , Ai,i 



1 1 

+ 



P{ij) + € P{ji) + e_ 
Dij^k = I dmidms . . . dmfc-i {irrii) {11111112) . . . {mk-2mk-i) {nik-ij) Dij,k , 



1=0 



n 



1 



n p{mri) +rej ^ P{msj) + {k - s)e^ 



mo = 1, ruk^ J. 



Up to /c = 2 we have 



hm {v,mT, -T)\v^) = {v,\v,) + —g 
1 — >oo ri 



1 1 

+ 



P{ij) + € P{ji) + 



+ ( ) / dm{im){mj) 



X 



1 



1 



1 



1 1 

+ 



_P{mi) + e P{ji) + 2e P(mj) + e P(ij) + 2e P(mj) + e P(mi) + e 
The formulas in A.) and B.) and C.) can also be written as 
hm U,{0,-T)\vj) = 

1 „ 1 



00 / ■\ k 
-I 



^ V ^ ; ^ i(^o - + ^e^'^^^Ki^o - + {k- 1)6 

1 „ 1 



int 



(i/0-£^°) + 26"'"*i(i/0-£^°)+6 



lim (^;,|C/,(T,0) = 



T-»oo 



{Vi 



CX) / .\ A; 
2 



k=l 



h 



int " 



:{Ho - E^) + {k - l)e ^{Ho - E^) + ke 
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lim {v,mT,-T)\v,) 

1 -^OO 



{Vi 



X 



oo k 



k=l 1=0 
1 



X 



-ffint • • • 



(6^ 



Below we present perturbation series up to second order for the components of the in 
and out states (JTOl) and (fTTl) . for the S-matrix elements (fT3il . for the energy eigenvalues 
(fT^ and ( JT5l) . for the phase factor ( JT2l) . and for the normalization constant (1201) . using 
the formulas written above. {vi\vj) are treated formally in the calculations (see 14. in 
section 1^751) . 6ij denotes {vi\vj) /{vi\vi). 

D.) For the in state (fTOl) we found 



\im^U,^,{-T){v,\UM-T)\v,) = {v,\v,)+g{^^]c,,^i,.+g^ 1^ ] c,,, 2,- + 0((?^), (69) 



where 



(70) 



dm 



{im){mj) 



[P{tj)+2e][P{mj) + e] 



+ / dm 



{jm){mj)5ij 
4e[e - P(mj)] 



E.) For the O'ui state (fTTl) we found 



dm 



{jm){mj)6ij 
4e[e + P(mj)] 

(71) 



^lim n,,,(T)*(t;,|f/,(T,0)|i;,) = (i^il^;,) + ^? [ ) + ( T ' ^^^'2,+ + 0(/), (72) 



where 



P(jz) + e ' 



(73) 



dm 



(im) (mj) 



[P(jz) +2e][P(mi) 



+ / dm 



{im)(jni)5ij 
4e[e - P(mi)] 



dm 



(zm) (mi) 6; 



{ii){ij) 



4e[e + P(mz)] 
(74) 
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F.) For the S-matrix elements ffT^ we obtained 



where 



lim n,,,(T)*n,-,(-T)(t;,|f/,(T, -T)\v,) 

T^oo 



dij,i = {ij){l-Sij) 



(75) 



(76) 



lij,2 



dm {im)(jrLj) 
1 



+ 



+ 



1-5, 



1 ^mi 



_P{ij) + 2e P(mj) + e P{ji) + 2e P{mi) + e P{mj) + e P{mi) + e_ 

{im){mi)6ij /nn^\/n^o\^. . 



+ / dm 



2e[e - P{mi)] 



2t[t + P{mi)] {vi\vi)2e'^ {vj\vj)2e'^ 



{v,\vi)e[P{ji) + e] {v,\v,)e[P{tj) + e] ' 
G.) We obtained the following formula for the eigenvalues of the in states: 

\imT^MHU,{0,-T)\vi) 



limT^oo(Wi|t^e(0, -T)\vi) 

(vi\vi) hj {vi\vi)[P{mi) + e] 



For the eigenvalues of the out states we obtained 

\iuiT^o.{v,\HU,{S),-T)\vi) 



E. 



i,out,e 



(9) 



\imT^^{vi\U,{0, -T)\vi) 

po I , (^0 I ,2-^ f {tm){rm){l-Sim) . 



H.) For the phase factor (fT^ we obtained 



lim n, , 



-T) = l + g- ^ ' 



+0{g'). 



h {vi\vi)e 
ii)"^ f {im)(jrLi) 



2e2 {vi\viY 



dm 



4e[e - P{mi)]{vi\vi) 



+ / dm 



[im) [mi) 



(77) 



(7J 



(79) 



4e[e + P(mz)](t'i|?^*) 



^0) 



22 



We found that limx^oo ^i,e{—T) = limT^Qo^i,e(T)* up to second order. 
I.) For the normahzation constant (l20l) we obtained 
lim ^JXi^^{-T) = {vi\vi) 

and hniT^oo \/ Xi^e{—T) = hniT^oo \/Xi^jT) up to second order. 



+ / dm 



{im){mi) 
4e[e - P{mi)] 



+ / dm 



{im)ijni) 
4e[e + P{mi)]_ 



(81) 



B Regular izat ion of the eigenvectors 

It was mentioned in section 12.31 that comphcations may arise when one calculates 0131) 
using perturbation theory or other methods, due to the fact that often {vi\vi) is not a 
finite number and has to be regarded as a distribution. One way to tackle this 

situation is to introduce a regularization of the eigenvectors \vi). In this appendix we 
outline such a regularization method. 
One takes superpositions 

|vi,/i)= / da f^{i,a)\va), (82) 



where /z is a regularization parameter and /^(i, a) are suitable functions of i and a with 
the property that 

{vi,iJ,\vj,ij) = I da f^{i,ayf^{j,a) (83) 



is finite for all i and j. In the limit /i — > the original \vi) vectors should be recovered. 
In particular, lim {v i, fi\vj, fi) = should hold, where S{i,j) is the Dirac-delta 

distribution. 

In the various formulas, in particular in f|T0|) - f|T5|) . one should replace and \vj) by 
\vi, fi) and \vj,fi). After this replacement the calculations can be carried out, keeping ^ 
and e finite. Then one should take the fi limit and, subsequently, the limit e ^ 0^. 

In the case of the S-matrix elements, for example, one has 



Si, = lim lim lim n,,,,^(T)*n,- ,,^(-T)(t;,, /i|f/,(T, -T)|t;,-, //), (84) 

6^0+ ^J,^0 T^oo 



where 



rrr^ V f^\Ue{T, \Vi, fl) {Vi, /i f/e 0, T \Vi, fi) 
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